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ABSTRACT 

In this paper we examine properties of the variable source Sgr A* in the near-infrared (NIR) using a 
very extensive Ks-band data set from NACO/VLT observations taken 2004 to 2009. We investigate the 
variability of Sgr A* with two different photometric methods and analyze its flux distribution. We find 
Sgr A* is continuously emitting and continuously variable in the near-infrared, with some variability 
occurring on timescales as long as weeks. The flux distribution can be described by a lognormal 
distribution at low intrinsic fluxes (< 5 mJy dereddened with Ak s = 2.5). The lognormal distribution 
has a median flux of wl.l mJy, but above 5 mJy the flux distribution is significantly flatter (high flux 
events are more common) than expected for the extrapolation of the lognormal distribution to high 
fluxes. We make a general identification of the low level emission above 5 mJy as flaring emission and 
of the low level emission as the quiescent state. We also report here the brightest Ks-band flare ever 
observed (from August 5th, 2008) which reached an intrinsic Ks-band flux of 27.5 mJy (ttiks = 13.5). 
This flare was a factor 27 increase over the median flux of Sgr A*, close to double the brightness of 
the star S2, and 40% brighter than the next brightest flare ever observed from Sgr A*. 
Subject headings: accretion, accretion disks — black hole physics — infrared: general — Galaxy: 
center 



1. INTRODUCTION 

The very center of our galaxy houses the variable 
source named Sg r A*, first discovered in the radio as a 
compact source (jBalick fc Brownl 119741 ). The fact that 
this source is motionless (to bette r than 1 km/s) at 
the d ynamical center of the galaxy (jReid fc Brunthalerl 
2004), and its coincidence with the common focus of el- 
liptical orbits of stars tracked over the last decade and 
a half, clearly ass ociates it with a supermassive black 
hole of 4 x 10 6 M^ ([Schod el et al.l[2002l: (Ghez et al.H 2008t 
IGillessen et al.ll2009h . 

Sgr A* has been detected across the electromag- 
netic spectrum, at radio, submm, NIR and X-ray wave- 
lengths. At NIR and X-ra y wavelengths dGenzel et al.l 
2003; IBaganoff et al1[200l the emission is highly vari- 
able (factors up to ftT60 and 27 in the X-ray and NIR 
respectively; iPorquet et al. I [20081 this work) compared 
to the comparatively steady emission at longer wave- 
lengths. NIR peaks are detected more often than in 
the X-ray (peaks occur «1 and 4 times a day for X-ray 
and NIR variable emission, respectively; Baganofi 120031 : 
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lEckart et~a l. 2006a). Some NIR flares ha ye been detected 
witho ut any accompanying X-ray flare (iHornstci n et al.l 
120071 ). However, when both NIR and X-ray exhibit 
increases in emiss ion, the peaks in emission occur si- 
multa neously (e.g. . lEckart et al.ll2004tlDodds-Eden et al.l 
[2009h . 

The near-infrared lightcurves from Sgr A* ex- 
hibit ~1 hour long increases in emission that are 
often called 'flares' in the literature. A num- 
ber of these have exhibited very suggestive sub- 
structural features with timescales of ~ 20 minutes 
dGenzel et al.l 120031: lEckart et all I2006bb iTrippe et al.1 
120071: lEckart et al.l I2008al lbl: IDodds-Ed eVeta!T T2009') , 
possibly quasi-periodic oscillations (QPOs). However the 
existence of QPOs and even the use of the term flare to 
describe t he NIR v a riabilit y of Sgr A* has been ques- 
tion ed bvlDo et all (l2009aD (see also IMever et al.l 120081 
and IMever et al.f l2009) who argue that there is no true 
quasi-periodicity, just a variability process with a fea- 
tureless red noise power spectrum (e.g. a power spec- 
trum P(f) ~ /~ 2 where / is frequency). A stochastic 
source with a red noise power spectrum has higher vari- 
ability at longer timescales and could potentially be re- 
sponsible for the structures on longer timescales seen in 
the real lightcurves. The authors suggest that apparent 
flare peaks may simply be the highest observed flux ex- 
cursions in such a purely stochastic source and are not 
isolated events. 

A main reason for the two rather contrasting interpre- 
tations of the variable emission from Sgr A* has been 
that the nature of the faint emission from Sgr A* and its 
relationship to the high flux emission is uncertain. The 
NIR emission from the Galactic Center is dominated by 
the central cluster of bright stars, and adaptive optics at 
8-meter class telescopes is required in order to separate 
the faint source Sgr A* from the closest S-stars (even at 
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this resolution Sgr A* is still on occasion confused with 
a relatively bright star). Additional, faint stars may be 
present very close to Sgr A* which have not yet been 
tracked and identifie d as stars from astro metric moni- 
toring programs (e.g. Gill essen et all 120091 ) . While the 
dramatic high flux variability can be unambiguously at- 
tributed to the black hole, when a faint source is detected 
at the position of Sgr A*, it is not necessarily clear that 
the source is Sgr A*, faint stars, or a combination of both. 
Accordingly, it is not clear whether Sgr A* continues to 
emit at all at low fluxes, whether it exhibits a 'quiescent 
state' (a non- or weakly active low state), or whether the 
low flux emission continues to vary constantly with the 
same statistical properties as the high flux emission. 

In addition to this, an unbiased overview of the prop- 
erties of the near-infrared emission from Sgr A* can 
be difficult to obtain from the published literature be- 
cause of publication bias (bright events have individ- 
ual interest and are often published alone). However, 
some studies have looke d at the statistical properties . 
lYusef-Zadeh et al.l (|2006l ) and lYusef-Zadeh et all ([2009f ) 
presented lightcurves and flux distributions for Sgr A* 
for about ^11 hours of 1.6/rai and ~32 hours of 1.7/im 
data observed with NICMOS on the Hubble Space Tele- 
scope (HST). With the resolution of the HST the close 
stellar sources are not as well separated from Sgr A* as 
with the VLT or Keck Telescopes, and the stars S17 and 
S2 overlap with the Sgr A* source in these observations. 
The flux distributions were fitted with a Gaussian at low 
fluxes, which was attributed to the observational noise on 
constant sources (the contribution from S2, S17 and pos- 
sible quiescent emission) and a power-law at high fluxes, 
which was attributed to transient flares. The best fit 
models implied that Sgr A* was active (above the noise 
a t low levels) more than 40% of the time. 

IDo et all (|2009aT ) presented an analysis of six nights of 
K'-band (and one L'-band) observations at the Keck Ob- 
servatory, using an unbiased set of observations taken be- 
tween 2005 and 2007. A source at the position of Sgr A* 
was always detected in this dataset, with an estimated 
maximum 35% contribution from stellar contamination. 
These authors reported that the source Sgr A* was con- 
tinuously variable, based on the larger variance of Sgr A* 
compared to stars of similar brightness on five of the six 
K'-band observation nights. This was the data set used 
to investigate timing properties of Sgr A* in which it 
was claimed the data set was consistent with a feature- 
less red noise power spectrum with no quasi-periodicity. 
However, with a sum duration for the K-band observa- 
tions of about 12.1 hours, this data set did not sample 
well the higher fluxes of Sgr A*, i.e. the source was 
relatively faint compared to publications where variable 
emission with s uggestive quasi-per i odic structure have 
been reported (iGenzel et al.1 120031 : iTrippe et all 120071 : 
lEckart et all l2008al for a comparison of Ks-band peak 
emission from the literature see T able [J in the Ap- 
pendix ) Although the studies of lYusef-Zadeh et all 
|2006l ). lYusef-Zadeh et al.l (pool ) and IDo et alT ipOOQa) 
have gone some way towards understanding the statisti- 
cal properties of Sgr A* in the near-infrared, there has 
not yet been a study on a very large, unbiased dataset of 
the variability of Sgr A* where the rare high fluxes are 
also well sampled. 

In this paper, we analyse the Ks-band flux distribution 



of Sgr A* for the years 2004-2009 from 117 observation 
nights carried out with the VLT in Paranal, Chile with 
the aim of seeking the flux-dependent characteristics of 
the variability of Sgr A* at both high and low fluxes. We 
do this through investigation of the flux distribution of 
Sgr A*. The dataset of this p a per is ^12 times larger 
than the data set of IDo et all ()2009af) . In order that 
our results might be compared with other publications, 
we give a summary in Table [T] (in the Appendix) of the 
brightest and/or notable Ks/K'-band variable emission 
reported in the literature, as well as the faintest values or 
upper limits. All are reported in the literature with dif- 
ferent calibrations/extinction values, and we have done 
our best to determine the corrections to scale them to 
the calibration and extinction values used in this paper. 

In Section [2] we present our observations and the re- 
sults of using two different methods of photometry: a 
six year lightcurve (aperture photometry), and the 2009 
data set (PSF photometry). In Section [3] we present the 
flux distribution of Sgr A* and our results from various 
model fits to the flux distribution. In Section [4] we dis- 
cuss our results in the context of two variability states for 
Sgr A*: a low-level lognormally varying quiescent state 
and sporadic high flux flares. We summarize in Section 

m 

2. DATA 

Since 2002 we have observed the Galactic Center 
with the near-infrared adaptive optics- as sisted diffrac- 
tion limited imager NAC O at the VLT ([Lenzen et all 
120031 : iRousset et al.ll2003l ) . Much of the observation time 
2002-2009 was spent concentrated on the central few 
arcesconds, measuring the positions of the S-stars and 
monitoring Sgr A*. By now we have amassed a large 
dataset with which we can investigate statistical proper- 
ties of Sgr A*'s variability. While we utilised L' (3.8/im), 
Ks (2.18//m) and H (1.66/im) bands, by far the largest 
proportion of our data is taken in Ks-band, collected in 
either 13 mas pixel scale imaging, 13 mas pixel scale po- 
larimetry, or 27 mas pixel scale imaging mode. 

The presence of many stars close to the central black 
hole complicates the attempt to acquire accurate pho- 
tometry for the near-infrared Sgr A* source. Sgr A* 
is usually fainter than the surrounding S-stars, and the 
stars move. A star can on occasion get so close to Sgr A* 
in projection that the two sources bec ome confused. For 
example, the ~16 ye ar orbit star S2 (jGhez et all 120081 : 
iGillessen et all [2009) was confused with Sgr A* as it 
passed pericenter in 2002. More recently, the mxs ~ 
15.8 mag star S17 was confused with Sgr A* in 2006- 
2008. 

We use two different methods in this paper: 

1. We first attempted to obtain as large a homoge- 
neous, unbiased, dataset as possible, in order to 
obtain a best overview of the statistical proper- 
ties of the variability of Sgr A*, and to address in 
particular the variability at (the rarer) high fluxes. 
To do this we used aperture photometry because 
this allowed us to address the most important com- 
plication to the photometry of Sgr A* in 2004 to 
2009 - the star S17, confused with Sgr A* in 2006- 
2008. To deal with S17 we used a two- aperture pho- 
tometry method to determine the combined flux of 
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TABLE 1 

Flux comparison of bright/notable/faint Ks-band fluxes from Sgr A* reported in the literature 



Reference 


Obs. 
Date 


Reported flux 
[dered mjy] 


Photometric cal. 
rcscaling factor 


Extinction 
rescaling factor 


Rescaled flux 
[dered mjy] 


Brightest States 




Peak fluxes 








Gen+03 


15 Jun 2003 


13.2 (10.5+2.7) 


1.10 


0.76 


11.0 


Gen+03 


16 Jun 2003 


10 (7.3+2.7) 


1.10 


0.76 


8.4 


Mey+07 


6 Oct 2003 


22 a 


1.2 


0.76 


20.1 


Tri+07/Mey+06 


31 May 2006 


16/23 (+S17) 


- 


- 


(16.7) 13. 5 b 


Hor+07 


31 Jul 2005 
2 May 2006 
17 Jul 2006 


11.6 

26.8 
6.8 


~1.06 c 


0.52 


6.4 
14.8 
3.7 


Do+08 


3 May 2006 

20 Jun 2006 

21 Jun 2006 

17 Jul 2006 

18 May 2007 
12 Aug 2007 


0.8 
0.65 
0.4 
0.3 
0.6 
0.2 


~1.06 c 


» 


8.5 
6.9 
4.2 
3.2 
6,1 
2.1 


Eck+08 


15 May 2007 


24 (+S17) 


0.76 


0.76 


(13.9) 10.7 


this work 


5 Aug 2008 


30.7 (+S17) 






(30.7) 27.5 


Faintest States 












Hor+02 


9 May 2001 


0.09 ± 0.005 mjy 
(upper limit, not dered) 


~1.06 c 




0.95 + 0.05 


Sch+02 


5 Jun 2006 


2+1 mjy 


? 


0.76 


1.5 + 0.8 


Do+08 


2006: May 3, 
Jun 20, 21, Jul 17 & 
2007: May 18, Aug 12 


0.192 mjy 
(median flux, not dered) 

0.082 + 0.017 
(faintest flux, not dered) 


1.06 d 
1.06 




2.0 
0.9 + 0.2 


Sab+10 


23 Sep 2004 
23 Sep 2004 


2.4 mjy 
(upper limit, 
no stellar contam.) 
0.9 mjy 
(upper limit, 
full stellar contam.) 


0.88 e x 0.90 f 


0.76 


1.4 
0.5 



Note. — Reported fluxes from the literature for high and low flux states of Sgr A* in the literature, rescaled to match the 
photometric calibration and extinction used in this paper (where S65 has rriKs — 13.7 and Aks — 2.5; see text for details). For the 
rescaled fluxes, brackets denote the raw observed flux (including S17) and the value without brackets the S17-subtracted estimate 
of the intrinsic flux. Note that without more detailed analysis we have only been able to reasonably account for and subtract the 
contribution of S17, and not any fainter stars: thus in the rescaled fluxes quoted here there may still be order of 1 — 2 mjy stellar 
contribution to the flux. 

a though the lightcurvc shows a peak fl ux more like ^24 m jy, this is likely an overestimate since the peak was only observed in one 
polarization filter (45°; sec Figure 2 in IMever et all (1200711 . We take instead 2 2 mj y, which appears to be a better estimate of the 
integrated flux (i.e. ^450 + -F1350), according to the modeling of Meyer ct al. (2007). 

b we quote the value from our own photometry; this observation night was part of our 200 4-2009 dataset. 
c Scaling for -0.06 offset in absolute photometric calibration, assuming same calibration as I Do et all fl2009b|) . 

d Scali ng fo r -0.06 offset in abso lute phot ometric calib r ation (as determined from comparing magnitudes reported in |Do et all 
d2009bl) and lGillessen et al.| 1)20091) , see also lSabha et all d2010D . 

e Sca ling for 0.14 m ag offset in absolute photometric calibration (as determined from comparing magnitudes reported in lSabha et al.| 
120101 and IGillessen et al.|[2009|) . 
^ Scaling for apparent different zcropoint. 
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Fig. 1. — Lightcurve of Sgr A* 2004 to 2009 from NACO Ks-band observations. The top panels shows the lightcurve of Sgr A* + S17, 
produced with aperture photometry, versus time with all gaps longer than 0.1 day removed. The second panel shows the lightcurve produced 
from the same data and with the same method for S7, a star of similar flux to the Sgr A*+S17 when faint. The source Sgr A*+S17 is more 
variable at low levels than the comparison source S7. In the lowest panel we show the detrended lightcurve, computed from subtraction of 
difference of the median value from each year, and the year of lowest median value (2006). We do this because there is good justification 
that the longest timescale trends of the source are dominated by additional faint passing stars (and in the case of 2004, extra flux in the 
aperture from the halo of S2 which was much closer in that year). Subtracted offsets are roughly 2.2 mjy in 2004, 0.9 mjy in 2005, mjy 
in 2006, mjy in 2007, 0.3 mjy in 2008 and 0.6 mjy in 2009. 
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Sgr A* and the star S17 in all years 2004 to 2009 
with comparable accuracy. This is important be- 
cause to understand the flux distribution we have 
to understand the photometric errors at low fluxes. 
If all years are measured with comparable accuracy 
the long-term errors can be well approximated with 
a Gaussian (with width dependent on flux; see Sec- 
tion [3TTl]). We did not include the Ks-band data 
from 2002-2003 because the star S2 was so close to 
Sgr A* during this time that it contributed to the 
flux measured with this method. Since S2 is much 
brighter than S17 the flux of Sgr A*+S17 could not 
be measured with an accuracy comparable to the 
other years. 

2. We secondly looked at a subset of the data in 
more detail in order to address the nature of the 
variability at low fluxes. Determining the na- 
ture/existence of low level variability (and in par- 
ticular to distinguish it from observational errors) 
requires more precise photometry than could be 
achieved with aperture photometry, in particular 
with the inclusion of SI 7. We address the ques- 
tion of the nature of the low level variability by 
analyzing the (high quality) 2009 data, in which 
S17 is not confused with Sgr A*, in greater detail 
with a PSF-fitting p hotometric method Starfmder, 
iDiolaiti et ail I2000D . The flux errors for this 
method are much smaller and allow us to distin- 
guish with more certainty between true variability 
and flux errors at low fluxes; this data set is not 
completely unbiased since Sgr A* was only reliably 
detected in 87% of the selected images. 

By using these two methods in combination, we can over- 
come the individual difficulties of the methods and piece 
together a consistent picture of the near-infrared vari- 
ability of Sgr A*. 

2.1. Aperture Photometry of Ks-band data 2004-2009 

Our 2004-2009 Ks-band data set consists of « 12000 
image s3 The data were reduced in the standard way by 
applying a sky subtraction, a flat fi eld correction and a 
hot/dead pixel correction (see e.g., iTrippe et al.i r2007^ . 

We found it necessary to apply a quality cut to elim- 
inate the worst data and ensure we obtained a homo- 
geneous data set on which we could perform accurate 
photometry. This quality cut was carried out by eye, 
and the criteria for elimination included: 

(i) the two calibration stars were not in the image (this 
happened rarely, and only in the polarimetric data) 

(ii) image ghosts close to Sgr A* or PSF artefacts 

(iii) simply a bad quality image mostly corresponding 
to low Strehl ratio data resulting from poor atmo- 
spheric conditions. The quality of images for this 
criterion was judged by the general visibility of S- 
stars, not by the visibility of Sgr A* to avoid bias 
towards bright fluxes from Sgr A*. 

9 proposal IDs: 072.B-0285, 073.B-0084, 073.B-0775, 073.B- 
0085, 271.B-5019, 075.B-0093, 077.B-0014, 078.B-0136, 077.B- 
0552, 078.B-0136, 082.B-0952 and Large Programs 179.B-0261/.B- 
0932 and 183.B-0100. 



We additionally eliminated 854 images from July in 2009 
which were taken in a triggered mode, so that our dataset 
remains unbiased and representative of the overall vari- 
ability of Sgr A*. The remaining data set totaled 6774 
images which were used to obtain a lightcurve for Sgr A*. 

To extract fluxes for Sgr A* as well as for several con- 
trol stars, we carried out aperture photometry on each 
image. However, if we used a standard circular aperture 
centred on Sgr A*, we would have a varying contribution 
from S17 (the brightest star to be confused with Sgr A* 
2004-2009) to the lightcurve as the star moves through 
the aperture during 2006-2008. A second star within 
the aperture, but not centered within the aperture also 
increases the error on our measurements of the flux of 
Sgr A*. We would obtain more accurate results if S17 
could somehow also be well-centered in the aperture at 
all times, as well as Sgr A*. 

To solve the problem of S17 we thus used a two- 
aperture method with two circular sub-apertures, one 
centered on Sgr A* and one centered on S17 and mea- 
sure their combined flux, averaging the results where the 
two sub-apertures are each 40, 53 and 66 mas in size. We 
used this method for all data 2004 to 2009, so that we al- 
ways measure the combined flux of S17 + Sgr A*. While 
there may be additional effects due to confusion with 
other, fainter stars, at first order our method ensures 
that we measure the fainter fluxes of Sgr A* with simi- 
lar measurement errors from year to year than we would 
have if we neglected to include a sub-aperture about S17. 

To further ensure the self-consistency of our dataset, 
we used the same set of calibrator stars for calibration 
of each image. This restricted us to only two suitable 
calibration stars, S30 and S65, due to the small field 
of view of some (especially polarimetric) images. To 
measure the raw counts of control stars and the cali- 
bration stars we used the same aperture as for Sgr A* 
and S17, but centered the star in just one of the two sub- 
apertures. For each image, the raw counts of Sgr A*/S17 
and the control stars are then divided by the counts of 
the two calibrator stars and the calibrated counts are 
then flux calibrated by scaling relative to the median of 
S65 for t he observation night, u sing mxs = 13.7 mag for 
S65 (e.g. iGillessen et al.l 120091 ) . Our photometric cali- 
bration is also consistent with lSchodel et ai1 [2010 (S65: 
rriKs = 13.64 ±0.02 ±0.06 mag). To convert magnitudes 
to flux es (mjy) we used the zeropoint from iTokunagal 
(2000) (ijiks = mag corresponds to 667 Jy). 

We computed the expected positions of Sgr A* and of 
the control sta r s usin g the orbital (polynomial) fits of 
IGillessen et al.l (|2009l ). For polarimetric data, we added 
the raw photon counts recorded in the apertures for all 
targets in both ordinary and extraordinary i mag es first 
in order to recover the unpolarized intensities^, before 
the flux calibration was carried out. We dereddened the 
final flux values using A Ks = 2.5 (ISchodel et al.l I20T0I 
Fritz 2010 in prep.). Note that the absolute photomet- 
ric calibration and extinction correction is not important 
for the analysis of relative fluxes presented in this paper, 
where we have focused on maximising the long-term rel- 

10 When using NACO in polarimetric mode, images in two po- 
larization filters are recorded at once on the detector by the use 
of a Wollaston prism plus mask: an ordinary image (in a chosen 
polarization filter) and an extraordinary image (at a polarization 
of 90° from the ordinary image). 
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Fig. 2.— The August 5th, 2008, lightcurve. Sgr A* + S17 is 
shown in blue, and the lightcurve of the star S7 is shown for com- 
parison. The flux has been dereddened with Ak s = 2.5. This flare 
is the brightest Ks-band flare that has ever been observed; the 
source reached an intrinsic peak flux of 27.5 mjy, a factor Ril.7 the 
flux of the star S2. While the high flux excursion is preceded by 
what looks like a flatter, 'background' level of emission, this is how- 
ever, at~7 mjy, several mjy brighter than the low flux levels from 
most other observation nights, and indicates there was increased 
source activity even at this time. The long-term trend for 2009 has 
not been subtracted from the lightcurve shown in this plot, and 
the last six datapoints did not make the data quality cut for the 
2004-2009 flux distribution. 

ative accuracy of our photometric measurements. How- 
ever, where we compare with previous publications, we 
additionally scale the reported fluxes to the photometric 
calibration (as best as possible) and extinction of this 
paper. 

The lightcurve we produced in this way for Sgr A* is 
shown in Figure [T] The data presented here is represen- 
tative of the overall variability of Sgr A* since it is not 
biased towards the presence of obvious 'flare events' and 
the only data selection carried out concerned data qual- 
ity and photometric consistency as explained above. It 
is by far the most extensive dataset for Sgr A* that has 
been published: the observations we present here consist 
effectively a ^184 hour long lightcurve (if one subtracts 
from the total time all gaps larger than about 20 min, i.e. 
gaps not due to sky observations). The total amount of 
exposure time (i.e. the time spent with the shutter open, 
excluding overheads) is smaller, and totals ~72 hours, 
though 184 hours is a better indication of the amount of 
continuous coverage. 

Figured] shows the 2004-2009 lightcurve of Sgr A* and 
comparison star S7. On top of the more rapid variabil- 
ity, there is a general trend for the lowest fluxes of the 
lightcurve to wander from year to year. It is at its low- 
est in 2006. In 2004 the extra flux appears to be due 
to a combination of a faint staiFI. a contribution from 
the halo of S2, which was much closer to Sgr A* during 
that year, and S19, which was confused with S17. The 
lightcurve also shows some increased flux in 2008 and 
2009; and as we will show in Section 12.1.11 there is at 

11 identified as S62 in Sabha ct al. (2010), though we disagree 
with this identification (S. Gillessen, private communication) 



least one additional (previously unidentified) faint star 
which must also contribute to the increased flux in those 
years. 

Since we can generally identify that stars are responsi- 
ble there is good justification for subtracting these long 
(i.e. year) timescale trends from the lightcurve. For want 
of a better method of determining the stellar contribution 
to the lightcurves we subtracted the difference between 
the median of the data from each year and the minimum 
median value (in 2006). This brings the offset in the 
lightcurve to roughly 3.5-4 mjy, much of which can be 
reasonably attributed to S17 (between 2.9 and 3.3 mjy 
in the lightcurve can be attributed to S17; see Section 
I2~2l) . 

It is apparent from Figure [1] that Sgr A* is much 
more variable than the comparison source S7. Some of 
the flux excursions are much more dramatic than oth- 
ers, with fluxes above 8 mjy (~5 mjy intrinsic upon 
subtraction of S17). Many of these have been previ- 
ously published - for example, the second brightest peak 
of the dataset, in May 20 06, was publi s hed p reviously 
in iMever et all (|2006D and iTrippe et all (|2007f l and has 
16.7 mjy (^13.5 mjy intrinsic; see Table [T] in the Ap- 
pendix) . 

On 5th August 2008, we saw a particularly extreme 
event which can be seen as the most dramatic flux ex- 
cursion in the lightcurve in Figure [1] This event is also 
shown in higher time resolution in Figure [2] On this par- 
ticular night the source at the position of Sgr A* bright- 
ened by a factor > 4 in a period of ^40 minutes, reach- 
ing a peak flux of 30.7 mjy (intrinsic flux of « 27.5 mjy 
minus SI 7, a factor increase over the upper limits on 
the lowest fluxes from Sgr A* of «27). An increase in 
Ks-band emission of this brightness is unparalleled in 
the literature for Sgr A*: the second brigh test that has 
been published is a flare from October 2003 (jMever et al.l 
2007) which reached an intrinsic flux of ~20.1 mjy (see 
Table [1] in the Appendix). 

2.1.1. Faint Stellar Confusion 2008-2009 

To create the 2004-2009 flux distribution of Sgr A* we 
subtracted a trend from the data that we suspected to be 
due to faint stars. However the trend for higher fluxes in 
2008-2009 is not explained by the presence of any known 
(previously identified) stars. 

Here, however, we are able to present evidence that 
there was indeed one to two faint stars confused with 
Sgr A* in the years 2007 to 2009. In Figure [3] we show 
that there are deviations in the source position of Sgr A* 
measured from imaging mosaics taken between 2007 and 
2009. The systematic deviations imply there was at least 
one, previously unidentified, faint star confused with 
Sgr A*. In fact in new, very high quality observations of 
early 2010, it is possible to identify two new faint stars 
very close to Sgr A*. One or both of these stars may 
have undergone a close pericenter passage. 

From Figure [3] it appears that one of the stars was 
present in 2007 also, while our offset for 2007 does not 
argue for any extra star. However as we discuss in Section 
14.21 there appears to be a relatively steady (on the order 
of ^0.5 mjy) stellar contribution to the flux of Sgr A* 
(including the years 2006-2007), which, together with the 
subtracted offsets, is probably the result of exactly such 
faint stars moving slowly in and out of confusion with 
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Fig. 3. — Measurements of the source position in declination for 
Sgr A* measured from (red) Ks-band and (blue) H-band imag- 
ing mosaics between 2007 and 2009. The clear deviations show 
that there was at least one contaminating star that may have been 
undergoing a close pericenter passage. Note that the measured 
deviations are dependent on the intrinsic fluxes of Sgr A*, which 
is variable, and the contaminating star. The zero position of the 
coordinate system is uncertain on the order of ~2-3 mas; for the 
purposes of this plot the zero position of declination was fixed to 
the position for which the brightest flux from Sgr A* was recorded 
(the 27.5mjy, mj^ s = 13.5 mag, event from August 5th, 2008). 

Sgr A*: probably one of the newly discovered stars does 
indeed already contribute to the flux in 2007 (with the 
additional increase in flux in 2009 due to the second star) . 

It is somewhat a mystery why these stars have not 
been identified previously. The stars responsible for the 
observed deviations in the source position have not yet 
been identified in data prior to 2007 when they should 
have been resolved separately. Perhaps they have un- 
til now always been confused with other stellar sources. 
The star close to Sgr A* in 2004 (identified as S62 in 
iSabha et al.ll2O10L though we disagree with this identifi- 
cation) may be a candidate. 

2.2. PSF photometry of Ks-band data from 2009 using 
Starfinder 

We additionally carried out PSF photometry on the 
Ks-band data from 20 09 (13 mas pixel sc ale imaging 
only) using Starfinder (|Diolaiti et al.l I2000D , which we 
use especially for investigating the question of whether 
Sgr A* is continuously variable at low fluxes (Sec- 
tion I3.2j) . We use an automated progra m for r unnin g 
Starfinder on many images developed bv iRankl {2007). 
The measured flux es were again calibra ted on the S-star 
S65 (mxs = 13.7; Gill essen et al.l [20091) and the results 
are shown in Figure |U We note that due to the pres- 
ence of the two previously unidentified faint stars (Sec- 
tion [2IlT|), we improve our chances of detecting a source 
in most images, while at the same time the faintness of 
this additional stellar contribution ensures our flux mea- 
surements have very small error, such that this data set 
is optimum for investigating the variability of Sgr A* at 
low fluxes. 

To produce the lightcurve we first required that a faint 
S-star (S21, of ~1.3 mJy) was detected in a given im- 
age as a basic data quality cut. From these images we 
then examined the detected positions of sources in the 
near neighbourhood (~100 mas) of the nominal posi- 
tion of Sgr A* (see Figure [4j and selected sources that 
were within 20 mas of the mean position of this detected 



source (by this method a source was detected in 87% of 
the good quality selected images). Due to this method 
of source selection, the dataset we present here is not 
completely unbiased, due to the 13% of images in which 
no source was detected, but the fraction in which Sgr A* 
was not detected is small enough that the data set still 
serves well to assess whether Sgr A* is in general more 
variable than stars of comparable flux. 

From the figure, Sgr A* certainly appears most of the 
time more variable than the comparison stars of similar 
flux (S21 and S60 in Figure @|. S60, in particular pro- 
vides a very good comparison since it is also close to S31, 
a star of similar brightness to S17; the two sources S60 
and S31 also have a very similar separation to Sgr A* and 
S17 during 2009 (see Figure gj. The flux distribut ion of 
S60 is fit by a Gaussian with n = 1.309 and a = 0.098 
which does not give any indication that the photometric 
accuracy is signficantly worse for sources of this separa- 
tion and flux ratio (the more isolated star S21 is fit by a 
Gaussian with /i = 1.333 and a = 0.090). 

In general the Starfinder and aperture photometry 
fluxes agree very well, with a constant offset of ~ 2.9 mJy. 
However, there is a systematic increase in the offset be- 
tween the measured fluxes from the aperture photome- 
try method and the Starfinder method towards the end of 
2009, when the source(s) confused with Sgr A* has(have) 
moved the furthest to the south (see Section [2.1.ip . This 
is likely explained by Starfinder's one-PSF fit to the dis- 
torted PSF of two sources (which has the effect of de- 
creasing the flux measured with the Starfinder method) ; 
the sources are even beginning to become separated at 
the end of the year (Figure 2]). 

This method also provides us with an estimate for the 
flux contribution of S17 to the 2004-2009 flux distribu- 
tion. The offset of ~ 2.9 mJy is slightly fainter than the 
flux of S17 (3.3 ±0.2 mJy when separated from Sgr A* in 
the 2009 data using Starfinder). The lower total fluxes 
obtained with the aperture photometry method is proba- 
bly a result of the higher stellar background surrounding 
Sgr A* (i.e. in the background aperture) compared to 
the more isolated flux calibration star and/or a side ef- 
fect of the double aperture method. 2.9 mJy is then a 
low estimate for the contribution of S17 to the combined 
Sgr A*+S17 2004-2009 lightcurve, while 3.3 mJy from 
Starfinder photometry is a high estimate. Therefore we 
estimate the contribution of the flux of S17 to our 2004- 
2009 lightcurve produced via aperture photometry to be 
in the range 2.9 — 3.3 mJy. 

3. RESULTS 

3.1. Flux Distribution 

For the flux distribution, we define the detection fre- 
quency in bin k as 



frcq 



del. A: 



(F k+1 - FfcjE^At, 



(1) 



where Aij and Fi are the exposure time and flux of the 
zth image, respectively, there are N images and Fk and 
Fk+i denote the bin edges. The detection frequency in 
each bin is divided by the bin width so that the area un- 
der the measured flux distribution is equal to 1. In the 
flux distribution each lightcurve data point was weighted 
by its exposure time so that the flux distribution is not 



DODDS-EDEN ET AL. 



16 
14 
12 

>; io 



-J 



31 10 



12 1 13 1 1 48 1 1 1 11 | 19 1 1 18 1 1 1 11 1 39 1 1 1 1 1 2 1 1 i I 1 17 
lOOmin 

Sgr A* (ap phot) 
including S17 



i 

o.i 



Sgr A* 
(Starfinder) 



. (?\ i 
'V. * ' 



0.5 1.0 
time [days, gaps removed] 



1.5 



_ 6 

1 4 

tr> 2 



■ 



a - - -o a g r - a- _ _^ ^. _ _g_ _ J_ _ _i 

starfinder - ap phot for SgrA* 



_ 6 

>-, 

? 4 



«3 £ 



S31 
'••S60 



S17 



S31 



S60 



SgrA* 



SgrA' 



S31 and S60 
Mar-Oct 2009 



S17 and Sgr A* 
Mar- Aug 2009 



-50 S17andSgrA* 
Sep-Oct 2009 



200 250 
A x [mas] 



-50 
A x [mas] 



-50 
A x [mas] 



Fig. 4. — Top: Lightcurves from aperture photometry (described Section |2. Il l and from Starfinder for 2009 Ks-band, 13 mas pixel scale 
imaging data. We have not removed here the long-term trend of 0.6 mjy that we did for the 2004-2009 lightcurve. Large gaps in time are 
removed so that the data can be seen in better time resolution; the real gaps in time between the individual datasets are shown at the top 
of the figure in days (the first dataset is 68 days from Jan 1, 2009). The data were selected on the basis of whether the star S21 (1.33 mjy) 
was detected or not, and then only those images (87%) in which a source was detected within 20 mas of the position of Sgr A*. Data from 
September-October 2009, where the faint star confused w ith S gr A* begins to be separated (see bottom panel) is shaded in gray: this data 
was not used for creating the flux distribution in Section 13.21 Second from top: Median flux differences between aperture photometry and 
Starfinder measurements. There is a systematic increase in the offset between the measured fluxes from the aperture photometry method 
and the Starfinder method towards the end of the year. Third from top: Lightcurves of comparison stars: S31 (2.97 ±0.16 mjy) and S60 
(1.30 ± 0.11 mjy). S31 and S60 are of similar separation and flux ratio to S17 and Sgr A* (when Sgr A* is faint). Bottom: Detected 
positions using the Starfinder algorithm on 2009 Ks-band, 13 mas pixel scale imaging data. S31 and S60 are two stars of similar separation 
and flux to S17 and Sgr A* when it is faint. Between March and A ugust both S17 and Sgr A* are well-detected on the images; from 
September onwards however the contaminating star (see Section l2.1.1l l begins to separate from Sgr A* and the photometry is unreliable. 



biased towards observations with shorter time samplings. 
We chose a logarithmic binning with bins spaced at in- 
tervals separated by a factor of 10 ' 05 . The errors for 
each bin are computed as the square root of the expo- 
sure times in the bin added in quadrature (with the same 
normalization as in Equation [T]) i.e. 



c-(freq dctifc ) 



^ =1 (AU(F k <Fj <F k+1 )y 
(F k+1 - F^l.AU 



(2) 



The flux distribution of Sgr A* for the years 2004-2009 
as derived from aperture photometry is shown in Figure 
[6] with various model fits (described in the next Section), 
as well as for a comparison star. The parameters and fit 
results for the model fits are given in Table [2J 

3.1.1. Accounting for observational errors 

The main difficulty with fitting the distribution of 
fluxes is the influence of observational errors, which are 
important enough at low fluxes to widen the peak of the 
distribution from its intrinsic shape. 
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Fig. 5. — Upper panel: Demonstration of the effect of observa- 
tional errors on the flux distribution. We show the flux distribu- 
tion from the Starfinder 2009 photometric dataset compared to the 
flux distribution determined via aperture photometry for the same 
(2009 subset) dataset. The fluxes obtained from aperture photom- 
etry method include the additional flux of S17. The size of the 
bins is scaled proportional to the typical error at the peak of each 
distribution. The blue dashed line shows the best fit lognormal 
model to the Starfinder photometry distribution, convolved with 
the error law for 2009 Starfinder photometry data a = 0.06F 0,6 
(Equation [4} , which are almost negligible (compare to the intrinsic 
flux distribution shown in gray). For the same data, but analyzed 
with aperture photometry, the observational errors are important 
and have the effect of slightly widening and flattening the peak 
of the distribution. The red dashed line shown is not a fit to 
the flux distribution produced from the aperture photometry data, 
but rather, it is the best fit model for the Starfinder photometry 
with a constant flux offset and convolved with the best fit error 
law dete rmi ned from stars for 2009 aperture photometry data (see 
footnotc ll2l) . Lower panel: Comparison stars with fluxes close to 
the mode the flux distributions shown above, demonstrating the 
much greater variance in the flux distribution for Sgr A* than that 
expected from observational noise in th e photometric data set ob- 
tained with Starfinder (see Section [32}. 

We characterized the observational errors by fitting 
Gaussian functions to the flux distributions of 22 nearby 
stars of fluxes between w0.6 and 25 mjy. Fitting a power 
law we obtain 

logio aobs(F) = (-0.76± 0.08) + (0.5 ± 0.1) log 10 F. (3) 

similar to that of IDo et all (|2009afl who found a typi- 
cal dependence of 0- o bs(F) ~ 0.2F 3 mjy for their data. 
The ~0.5 power-law dependen ce of the errors i s consis- 
tent with a photon noise origin (IFritz et alJI2OT0h . We do 
not find any flattening of the power law dependence of 
the noise towards low fluxes with statistical significance. 

The 13 mas pixel data from 2009 is of much higher 
quality and for the Starfinder photometry has noise prop- 
erties 

logio <Jobs(F) = (-1.2 ± 0.1) + (0.6 ± 0.2) log 10 F. (4) 

The noise dependence was determined in the same way 
as for the 2004-2009 dataset through comparison of 22 



stars of fluxes 0.5 to 25 mJvP^l 

We include the effect of observational errors in our 
model fits by convolving the intrinsic flux distribution 
with a Gaussian of width a b s (F) , which has the depen- 
dence on flux of Equation [3] 



F 



0+crr 



(F) = J Po(F') 



1 



27TCT obs (F) 



exp 



-{F-F'f 
2a ohs (F) 2 



dF'. 



(5) 

where Pq is the intrinsic flux distribution (without ob- 
servational errors). The effect is only important at low 
fluxes < 6 mjy (see Figure [6]) and is responsible for 
broadening the low flux peak in the observed flux dis- 
tribution for Sgr A*. 

3.1.2. Lognormal model 

A lognormal fit to the 2004-2009 flux distribution from 
aperture photometry is shown in Figure |B1 with the best 
fit parameters given in Table [2] This model has intrinsic 
probability distribution Pq = P\ ogn where 



flogn(F) 



2na*(F-F b ) 



exp 



Qn(F ~ Fb) - I**) 

2(7 , 2 



( 6 ) 

where F represents the observed flux (measured in mjy 
here and elsewhere in this paper), and F > F b with F b 
a flux offset due to some constant contribution (e.g. S17 
+ contaminating stars). The parameters /i* and er* of 
the lognormal distribution have a natural analogy to the 
normal distribution when exponentiated: the source can 
be thought of as having a median flux of exp(/z*) with 
a multiplicative standard deviation of exp(<7„) (i.e. the 
interval exp(/z*)/ exp(<7*) to exp(/i») x (exper*) contains 
«68% of the probability). 

3.1.3. Lognormal model with flatter high flux tail 

The residuals at low fluxes of the lognormal model are 
probably just due to the low flux wing of the error dis- 
tribution which is not perfectly Gaussian (see Star S7 
in Figure [6]). The residuals for the lognormal model at 
high fluxes, on the other hand, can not arise from mea- 
surement error, and we obtain a much better fit to the 
flux distribution with a lognormal model which makes a 
transition to a flatter tail at high fluxes. The intrinsic 
probability distribution in this case is 



-Plogn+taiK-F 1 ) 



kP iog n(F) 

kP loga (F t )(^y 



F < F b 
F > F b 



F, 
F t 

where P\ gn(F) is the lognormal distribution of Equation 
[6j F t is the flux at which the distribution makes the tran- 
sition to the flatter tail, s is the power-law slope of the 

12 the increase in accuracy appears to be due to both the higher 
quality of 2009 data and the higher accuracy of the Starfinder 
method: for example, for the same data of 2009 used for the 
Starfinder data set, but analyzed with the aperture photometry 
method, the noise properties are 

log 10 <T(F) = (-1.0±0.1) + (0.5±0.2)log 10 F, 

which is a significant improvement on the overall 2004-2009 aper- 
ture photometry accuracy, but still not as good as the accuracy 
obtained with Starfinder for the same dataset. 
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Fig. 6. — Left and center: Lognormal and lognormal+tail model fits to the 2004-2009 flux distribution (see Sections 13. 1.2l and 13. 1.3( l. The 
lognormal+tail model with a transition to a flatter tail towards high fluxes is a significant improvement to the fit, at > 6cr confidence. Note 
that although the transition flux appears to be at about 8 mjy (since the contribution of S17 is always added to Sgr A* in this photometric 
dataset), the intrinsic transition flux is actually Ft = 4.6 ± 0.5 (see Tabled- The blue lines show the full model flux distribution including 
observational errors (o~(F) = 0.17 F 0,5 ; Equation 13) . while the dashed gray line shows the intrinsic flux distribution (see Section l3.1,l| l. 
Essentially all measurements below ~3.5 mjy occur due to observational errors. The best fit parameters for these models can be found in 
Table [2] Right: the flux distribution of comparison star S7, together with the best fit Gaussian. The Gaussian has fi = 4.9 and a = 0.48 
for a x^/doi = 37.2/8. The largest contribution to \ 2 ( on the order of 20 to 30) for both S7 and Sgr A* best fit models comes from the 
low flux wing of the distribution, where a Gaussian fails to be a good model for the observational errors. 



TABLE 2 

Fits to the observed flux distribution 



logn+err (see Section 13 . 1 .2D 



ln(mjy) 



ln(mjy) 



F b 
mjy 



X 2 /dof 



logn+tail+err (see Section 13.1.3ft 



ln(mjy) 



ln(mjy) 



mjy 



Ft 
mjy 



X 2 /dof 



Aperture Photometry on 2004-2009 data 

-0.12 + 0.04 0.92 + 0.03 3.75 + 0.04 79.7/18 



Aperture Photometry on 2004-2009 data 

0.05 + 0.06 0.75 + 0.05 3.59 + 0.06 4.6+ 0.5 2.70 + 0.14 47.8/16 



Starfinder Photometry on 2009 subset 
0.04 + 0.06 0.77+ 0.05 1.34+ 0.05 19.0/14 



Starfinder Photometry on 2009 subset 

0.07 + 0.06 0.73 + 0.06 1.31 + 0.06 4.8+ 0.9 2.1 + 0.6 15.3/12 



Note. — Best fit parameters and formal fit errors for the large (aperture photometry 2004-2009) and high quality (Starfinder 2009) 
photometric datasets. For description of the models and the parameters see Sections 13.1.11 13.1.21 and 13.1.31 Note that F t refers to 
intrinsic flux, while on the observed flux distribution the break appears to occur at a flux of Ft, + Ft. 



tail, and k is a renormalizing factoid such that the total 
probability is 100%. 

The flattening of the distribution towards high fluxes 
in the 2004-2009 aperture photometry flux distribution 
is significant at > 6a. However, the fits are not formally 
'good fits'. This is to the most part due to residuals in 
the lowest flux bins, most likely caused by the deviation 
of the observational errors from the plain Gaussian de- 
scription used in our model. If we instead only fit only 

13 k = 1/(| + |erf((lnF t - ^)/y/2b~^) + F t P logn (F t ) / (s - 1)) 



to fluxes > 3.5mJy (ignoring the lowest three bins), we 
obtain a more reasonable % 2 /dof for the lognormal plus 
tail model of 17.4/13 (while the best fit lognormal model 
for the same has 37.6/15). In this case the flatter tail of 
the distribution is still significant at > 4c. 

Interestingly, in the 2009 Starfinder photometric 
dataset the tail is also present with similar best fit pa- 
rameters, although in this smaller dataset (which only 
reaches fluxes of ~10 mjy) the tail has a significance of 
only 1.4<7. 
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3.2. Is Sgr A* a continuously variable source? 

What is not clear in the 2004-2009 data set is whether 
it is appropriate to fit only continuously emitting, contin- 
uously variable models to the flux distribution. A plau- 
sible alternative is a model in which the source of near- 
infrared emission emits only sporadically (i.e. exhibits 
'on' and 'off' periods). If this is the case observational 
errors alone (on the measurement of a constant compo- 
nent to the flux) should be responsible for observed vari- 
ations in the flux level during 'off' periods. The goal of 
this section is to determine whether there is too much 
variation at low levels for the data to be explained by 
observational errors in 'off periods' of a sporadic model. 

We use the following prescription for a sporadic model 
for the flux distribution: 

P +err(F) = PvaI J P (F' )P gauBS (F - F' - F b )dF' 

+ (l-p 

var ) -Pgauss 

(F - F b ) (8) 

where p vaT is the fraction of time the source spends in an 
'on' state. 

Using Pq(F) = P\ ogn (F) (Equation E]) and fitting to 
both the 2004-2009 aperture photometry flux distribu- 
tion and to the 2009 Starfinder distribution with the mea- 
sured noise laws of Equations [3] and |U we find in both 
cases the best fits to the distributions require p val > 92% 
and pvar > 96% (at 3er) respectively, i.e. the measured 
noise laws can not account for the observed amount of 
variability at low fluxes. 

To look at the question from a different angle, we inves- 
tigate what level of observational noise would be required 
to produce the observed variation at low fluxes. To do 
this we allow the normalization of the noise law to be a 
free parameter in the model fits, i.e. 

log 10 a obs (F) =A + 0.5\og 10 F. (9) 

We find A = -0.66 ± 0.02 and A = -0.79 ± 0.04 for 
the 2004-2009 aperture photometric datset and the 2009 
Starfinder dataset respectively. We also refit the noise 
laws to the stars (see Equations [3] and HJ) fixing the slope 
of the error law to 0.5 and obtain A sta rs = —0.84 ± 0.05 
and A^rs = -1.12±0.07 for the 2004-2009 aperture pho- 
tometric datset and the 2009 Starfinder dataset respec- 
tively. The best fit noise law to the 2004-2009 flux dis- 
tribution from aperture photometry allows p vaI — 68 ± 7 
but the normalization of the noise law is inconsistent with 
that of the stars at a level of 3. 3c. For the 2009 flux dis- 
tribution from Starfinder photometry the required noise 
law is even more unlikely, and is ruled out at a level of 
4.5<7. We conclude that the observed behavior of Sgr A* 
at low fluxes is best consistent with a continuously emit- 
ting, continuously varying source. 

3.3. Variations on timescales of weeks to months 

We also detect longer timescale trends in the low flux 
level of Sgr A* in the 2009 data (Figure [4} , such as an 
~0.8 mjy increase in the minimum level from the fourth 
to the fifth data sections, with a gap in time of only 12 
days, and a similar size decrease from the second to the 
fourth data sections with a gap in time of only 11 days. 

It might be that a passing star (such as the two new 
stars discovered in 2010, see Section [2. 1.1[) is responsible 
for these trends, either passing in and out of confusion 



very quickly, or even just moving within the Sgr A* PSF 
and distorting the flux: the further apart the two sources 
the lower the fitted flux. However these stars cannot 
move quickly enough for this: the change in flux is so 
large that in either case the star must move by at least 
^50 mas within a period of 12 days to distort the flux 
as much as observed, which at the 8kpc distance of the 
GC requires an average speed of at least 0.1c over a two 
week timescale. Even assuming the most extreme case 
- a star moving only in the plane of the sky and falling 
towards Sgr A* in a parabolic orbit at the escape speed - 
a star should take at least ^110 days (0.02c in average), 
to cover such a distance. A larger centroid shift between 
the datasets would also be expected if it were the case 
that a star were distorting the source PSF (~20 mas 
instead of the observed ^5 — 8 mas) . Furthermore, even 
if this could explain the trend between the fourth and 
fifth datasets, the same explanation would have to be 
invoked to again explain the variation between the second 
and fourth datasets, also on a timescale of only ^11 days. 

We therefore conclude that the observed trends on 
timescales of a few weeks in the lightcurve are most likely 
intrinsic. Though they can not be due to direct flux 
contamination from passing stars, the trends are possi- 
bly still related to the feeding of the accretion flow from 
stellar winds of the passing stars, which may trigger in- 
creased act ivity as the star passes orbital pericenter (e.g. 
lLoebl[200l . 

4. DISCUSSION/INTERPRETATION 

4.1. Two states of Sgr A* in the near-infrared 

We have shown in Section [3] that the flux distribution 
of Sgr A* is best described by a continuously emitting, 
continuously variable component with a lognormal flux 
distribution at low fluxes, but additionally that there is a 
significantly flatter tail above about 5 mjy, implying that 
there is something different about the high flux emission 
from Sgr A *. Put together, it seems justified to identify 

• emission below 5 mjy as quiescent state emission 

• emission above 5 mjy as flaring emission. 

This is the (phenomenological) definition by which we 
hereafter refer to either 'flares' or 'quiescent' emission in 
Sgr A*. We present how the two states may fit into a 
possible physical picture below. 

4.1.1. Quiescent State Emission 

That the low fluxes of Sgr A* are fit by a continuously 
present, variable component suggests that there is truly a 
quiescent stat^\ in Sgr A*, i.e. a low- level, lognormally 
varying component (factor ^2 in flux of the lognormal 
component corresponds to la variability). Our finding 
of continuously emitting, co ntinuously variab le compo- 
nent confirms the findings of lDo et"aLl (|2009a|) who first 
reported continuous variability of Sgr A* from the larger 
variance of the flux of Sgr A*, as compared to stellar 
sources, in 5 out of 6 K'-band observation nights. 

Lognormal distributions of fluxes are seen in X-ray 
binaries, where they are interpreted as a sign that 

14 a quiescent yet variable state, in a similar sense to the quies- 
cent state of longer wavelengths which nevertheless exhibits vari- 
ability (factor 20-100% in radio to submm regime). 
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multiplicative pro cesses are behind the variability (e.g. 
lUttlev et al.l 12005). A good example of a multiplicative 
process is the model of lLvubarskiil f|1997T ) where variabil- 
ity is produced through the inwards propagation of accre- 
tion rate fluctuations at different radii. These variations, 
produced in regions of the flow which themselves do not 
emit near-infrared emission, all combine multiplicatively 
as they propagate inwards. The end result of combined 
fluctuations is then imprinted on the variability pro- 
duced by the (say, near-infrared) emitting region of the 
flow. A very similar idea may work with magnetic tur- 
bulence, which has been shown in GRMHD/MHD simu- 
lations to be capab le of producing the factor 40-50% mil- 
limeter variability (iDexter et al J 1 20091 iChan et al.l 120091: 
IGoldston et al.H2005D . 

From where exactly the quiescent emission arises is 
however a difficult question. Nonthermal electrons are re- 
quired to emit in the near-infrared, and it's not clear how 
these should be related to the thermal electrons emit- 
ting at submm wavelengths, or how they should neces- 
sarily be arranged within the accretion inflow/outflow. 
lYuan et al.l (|2003l see also lOzel et al.l 120001 ) require a 
hybrid thermal/non-thermal electron distribution to ex- 
plain the spectrum of Sgr A* (in particular the excess 
radio emission at low frequencies). The models pre- 
dict that the same non-thermal electrons (with ~1.5% 
of the thermal energy) could produce emission in the 
near-infrared. Thus it might be interesting to compare 
the flux distribution of the NIR quiescent state emission 
(F < 5 mjy) with that of the low-frequency radio emis- 
sion, as well as correlations in timing properties, such 
as whether the medium-timescale trends (~2 weeks) ob- 
served in the 2009 data (Figure|4]) correlate with any sim- 
ilar long-term trends in the radio regime. More detailed 
comparison of the low-level NIR emission from Sgr A* 
with other wavelengths (e.g. submm) could also shed 
light on how non-thermal electrons are related to the 
thermal population. 

4.1.2. Flaring Emission 

It seems unlikely that the same variability process is 
responsible for both high and low flux emission from 
Sgr A*. There is no obvious reason for the underly- 
ing process, if multiplicative at low fluxes, to deviate 
at high fluxes in order to explain the flatter tail above 
4.6 ± 0.5 mjy. The flatter tail of the flux distribution in 
Sgr A* requires at the very least some transition within 
the source (to a flaring state), triggered around an in- 
trinsic flux of 5 mjy, such that either large fluctuations 
or their emission undergo some kind of runaway amplifi- 
cation. 

On the other hand the physical mechanism behind the 
flares could be unrelated to (though possibly still trig- 
gered by) the mechanism producing the low level vari- 
ability. The flare tail could arise from additional sporadic 
flare events which dominate the distribution of fluxes 
above 5 mjy. In this case 5 mjy is not a physically in- 
teresting flux outside the fact that it is the flux at which 
one is equally likely to observe either a flare event or a 
low level lognormal variation. 

Spontaneous magnetic reconnection is a good candi- 
date for a physical process that could give rise to spo- 
radic flares in Sgr A* (e.g. lYuan et ail 120031 : iDing et all 



2010). The decrease in magnetic field strength that 
would accompany the conversion of magnetic energy to 
high energy electron acceleration in a magnetic recon- 
nection event can explain the non-one-to-one rise in the 
lightcur ves of the so far brightes t NIR/X-ray flare from 
Sgr A* (|Dodds-Eden et al.ll2010l ). 

A number of other properties also indicate that there 
are two types of NIR emission from Sgr A*, with differ- 
ences between high and low flux emission, which include: 

• Spectral Index. There are indications that the spec- 
tral in dex at low fluxes i s redder than at higher 
fluxes dGhez et al.l 120041: lEisenhauer et al.l 120051: 
iGfllessen et al.l 120061: lYusef-Zadeh et al.l 120091 ). al- 
though this has been disputed by other authors 
IHornstein et~aTl (|2007D . 

• Polarization. Changes in polarizatio n in the ris- 
ing/ de caying phases of flares are seen (|Eckart et al.l 
2006b), and low level emission appears to be gen- 
erally more highly polarized ( ~ 30 — 40%) than 
high flux emission (~ 10%) (|Mever et all 120061: 
iTrirjpe et al.ll2007l ). 

• Association with X-ray flares. The best fit model 
(Lognormal+tail+ errors) to the 2004-2009 flux dis- 
tribution has 4.5% of the probability above the best 
fit transition flux of F t = 4.6 mjy, which, interest- 
ingly, is close to the X-ray flaring rate of ~4 - 7% 
(for flares occu rring once a d ay with a duration 
of 60-100 min; iBaganoffl 12001 . Estimates of the 
NIR flaring rate including emission below 5 mjy, 
on the other hand, typically put the NIR flaring 
rate much higher at » 40% (jEckart et all I2006al: 
lYusef-Zadeh et all 120061 ) . Perhaps the additional 
low-level component is even responsible for differ- 
ences in the NIR/X-ray peak flux ratios. 

There are a number of high flux Ks/K'-band events 
(> 5 mjy), not part of the 2004-2009 dataset pre- 
sented here but published elsewhere in the literature, 
which are also of >5 mjy intrinsic flux and can be 
identified as belonging to the flare tail of the flux dis- 
tribution. To list some of the brightest (see Table Q] 

in the Appen dix for details ): ~11 and 8.4 mjy (15, 

16 June 2003: IGenzel et all 12001 : 20.1 mjy (6 Octo- 
ber 2003: IMever et all 120071): 13 .5 mjy (31 May 2006; 
iTrippe et all 120071 iMever et all 120071 and part of our 
2004-2 009 dataset): 10. 7 mJv (15 Mav 2007:lEckart et all 
I2008af ): and 14.8 mjy (|Hornstein et al.ll2007l ). ~ 

4.2. Consistency with previous measurements of Sgr A * 
at low Ks-band fluxes 

Our results are reasonably consistent with the me- 
dian flux of 2.0 mjy over s ix separat e observ ation nights 
from 2006-2007 found by IDo et all (|2009aD (see Table 
Q] in the Appendix for the scaling factors applied to 
scale the reported values to our calibration/extinction). 
We find a median flux of « 1.6 mjy for our 2004-2009 
lightcurve with the long-timescale trend and S17 (as- 
suming 3.1 mjy) removed (note however that this still 
includes some constant contribution: the median flux of 
the lognormal component alone is 1.1 mjy). 

A number of previous measurements have put up- 
per limits on any long-term steady quiescent state at 
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near-infrared wavelengths which have been on the ~1- 
1.5 mJy level (iHornstein et al.ll2003 iSchodel et~aLll2007l: 
ISabha et al.ll2010l see Table [T1 in the Appendix). These 
measurements are consistent with the picture we have 
put forward here for Sgr A* at low levels, of a con- 
tinuously emitting, continuously varying source. For a 
lognormally varying quiescent state with the properties 
of our best fit model to the 2004-2009 flux distribution, 
an upper limit of 1-1.4 mJy in selected observations is 
not surprising: the source is expected to be at <1 mJy 
roughly 50% of the time. 

However, more likely the upper limits of 1-1.4 mJy 
includ e a stellar contribution to the flux. ISabha et al.l 
(|2010fl estimated an upper limit on the intrinsic flux 
from Sgr A* of F <0.5 mJy in the presence of a likely 
~0.8-l mJy stellar contribution. In the context of our 
best fit lognormal model, this low flux measurement is 
a rarer, but not unlikely occurrence, with a probabil- 
ity of - 1 6%E3 The . 8-1.0 mJy diffuse component re- 
ported by Sa bha et~aLl (|2010l) might be associated with 
the ~0.3-0.7 mJy offset we find in our fits to the flux 
distribution of Sgr A*, between the fitted constant com- 
ponent (3.6 m Jy) and the expe cted contribution of S17 
(2.9-3.3 mJy). IDo et al.1 (|2009a|) also estimated from the 
color of Sgr A* at low fluxes, that they most likely had 
^35% contribution from stellar contamination (0.7 mJy) 
for their observations in 2006. These results all seem to 
point towards an additional stellar component to the flux 
of Sgr A* with a constant flux of ~0.5 mJy (our data), 
though it does not necessarily always consist of the same 
stars from year to year. 

4.3. Comparison with X-ray binary variability 

Lognormal distributions of fluxes are also seen in 
other astrophysical objects, for example X-ray bi- 
naries (prototype Cyg X-l in the low hard state, 
lUttlev et al.l I2005T ). implying a connect ion between 
Sgr A * other low-accret i on ra te black holes Falck e et al.1 
(|2004r i: iMcHardv etaLl (|2006f ). The physical picture of 
the low-hard state of X-ray binaries is that - compared 
to higher luminosity flows - the inner optically thick ac- 
creti on disc disappear s and is replaced by a hot inner 
flow (|D one et al.ll2007l ). Sim ilar pictures, w ith hot inner 
flows, exist for Sgr A* (e.g. IQuataertj|200l . Sgr A* is, 
however, at the extreme low end of the accretion rate 
range of observable sources (the 'quiescent' state). Con- 
nections between Sgr A* and low-hard state X-ray bina- 
ries of higher accretion rate such as Cygnus X-l can shed 
light on whether or not the quiescent state is indeed a 
distinct state or whether there is simply a smooth con- 
tinuat ion of the har d state down to such low accretion 
rates (|Markofl1l2010l) . 

If timescale scales with mass, then the flares in Sgr A* 
of ^10-100 minute duration, correspond to ~ millisecond 
timescales in a stellar mass black hole binary. It is also 
true that millisecond-timescale flares are seen in Cyg X- 

15 We note that for the same dataset (September 23, 2004) as 
was used to compute the upper limit of 0.5 mJy, our aperture pho- 
tometry method yields a flux of 1.0 mJy with aperture photome- 
try (with S17 and the empirical background fi t subt racted). This 
compares favourably with that of Sabha et al. (2010, 0.8-1.0 mJy, 
scaled to our calibration) obtained with aperture photometry after 
PSF-extraction of all known sources close to Sgr A* including the 
star identified in that work as S62. 



1, though it was shown by lUttlev et al.l (|2005l) that the 
majority of these (and all those in the low/hard state) 
were not more frequent than expected in the context of 
the extension of the low-level lognormal distribution in 
this source to higher fluxes. There was one flare however, 
a 12.4cr event that occurred (in fact when Cyg X-l was 
in the high/soft state) that could not be accounted for in 
the context of the underlying variability and was posited 
by the authors to have a different physical origin. In 
this vein, it could be that the flares in Sgr A* are due 
to the source's "flickering attempts at o utburst a ctivit y 
from out of quiescence" , as sugg ested bv lMarkofl (pOlOl) . 

A connection between the near-infrared variability of 
Sgr A* and the variability of X-ray binaries has also been 
made previously from a long timescale break in the PSD 
of Sgr A* which can be associated with a break timescale 
in the low-hard state of X-ray binaries (Meyer et al. 
2009), though the single dataset used to constrain the 
high frequency power spectrum in this case may not be 
representative. The lognormal distribution for low-level 
fluxes makes however an independent case for an anal- 
ogy between the variability of Sgr A* and low-hard state 
X-ray binaries. 

4.4. Wavelength Dependence 

The flux distribu tions of lYusef-Zadeh et al.l (|2006l ) and 
lYusef-Zadeh et all (|2009l ) at 1.6/zm and 1.7/zm display a 
qualitatively similar shape to our Ks-band distribution, 
although these datasets are still much smaller than ours, 
and the low-level activity appears to be dominated by 
the observational noise. Nevertheless, with more data at 
1.6 and 1.7 fim. and modeling of the observational errors, 
it might be possible to obtain constraints on the colours 
of the low level and high flux components. In future 
work we intend to use the same techniques as we have 
used in this paper at Ks-band, in other bands (L' and 
H) to investigate in detail the dependence of the flux 
distribution with wavelength. 

4.5. Timing Analysis 

Our work is not the first to point out that the flux 
distri bution of Sgr A* is not w el l-describe d by a G aus- 
sian: lYusef-Zadeh et ail (|2006l ). IDo et al.l (|2009aD and 
Yusef-Zad eh et alT()2009f ) have all noted that the flux dis- 
tribution has a high flux tail. However, t he timing analy- 
ses which have been carr ied out to date (|Do et al.l f2009a: 
IMever et aL| [2008. 2009) in which Sgr A* was compared 
to simulated red noise lightcurves have nevertheless used 
simulated lightcurves which have had, by construction, a 
Gaussian flux distribution. Sgr A*'s intrinsic mean and 
variance are such that if one were to plot the simulated 
lightcurves scaled to Sgr A*'s mean and variance (for the 
2004-2009 dataset fj, = 1.3 mJy subtracting 3.8 mJy for 
S17 plus other faint stars, and a — 1.9 mJy), the model 
red noise lightcurves would have negative intrinsic fluxes 
muc h of the time whic h is rather unphysical. 

In IDo et al.l (|2009aD the high flux tail was excluded 
when fitting the flux distribution, and only fluxes 
<3.2 mJy were fitted with a Gaussian (see Table Q] in 
the Appendix for scaling factors; the reported value 
m IDo et all I2009al was a reddened flux of 0.3 mJy). 
Thus their model lightcurves best resemble the low level 
emission of Sgr A*. Indeed the entire dataset used by 



14 



DODDS-EDEN ET AL. 



IDo et all (|2009al) does not include much bright emission 
> 5 mJy, and thus their results really apply to the low 
level emission of Sgr A*. That the low level emission 
from Sgr A* is consistent with red noise fits together 
with our finding of the lognormal distribution of fluxes 
at low fluxes. Both of these aspects of t he variability 
are se en in X-ray binaries such as Cyg X-l ()Uttlev et al.l 
2005). However, given that the flux distribution is not 
Gaussian, it remains questionable that the power spec- 
trum has yet been reliably constrained from the Monte 
Carlo simulations of previous timing studies. 

In particular, high and low flux states of Sgr A* may 
require separa te timing analysis. In this respect it is in- 
teresting that I Mever et al.l (2008) finds a quasi-periodic 
signal with false alarm probability of only 2 x 10~ 5 (4.2c 
significance in Gaussian equivalent terms), if one consid- 
ers just the window 385-445 minute subset from 30-31 
July 2005, though the signal was deemed insignificant 
when the entire lightcurve from the night was analyzed. 

5. CONCLUSIONS 

We have determined the distribution of fluxes for the 
variable source Sgr A* in the near-infrared and thereby 
investigated flux-dependent properties of the emission. 
We summarize our main results as follows: 

• Sgr A* is continuously emitting and variable in the 
near-infrared. 
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